Free Isotropic-Nematic Interfaces in Fluids of Charged Platelike Colloids 
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Bulk properties and tree interlaces ol mixtures ol charged platelike colloids and salt are studied 
within density-lunctional theory. The particles are modeled by hard cuboids with their edges con- 
strained to be parallel to the Cartesian axes corresponding to the Zwanzig model. The charges ol the 
particles are concentrated in their center. The density lunctional is derived by functional integration 
ol an extension ol the Debye-Hiickel pair distribution lunction with respect to the interaction po- 
tential. For sufficiently small macroion charges, the bulk phase diagrams exhibit one isotropic and 
one nematic phase separated by a first-order phase transition. With increasing platelet charge, the 
isotropic and nematic binodals are shifted to higher densities. The Donnan potential between the 
coexisting isotropic and nematic phases is interred from bulk structure calculations. Non-monotonic 
density and nematic order parameter profiles are found at a free interlace interpolating between the 
coexisting isotropic and nematic bulk phases. Moreover, electrically charged layers form at the free 
interlace leading to monotonically varying electrostatic potential profiles. Both the widths ol the 
free interlaces and the bulk correlation lengths are approximately given by the Debye length. For 
fixed salt density, the interlacial tension decreases upon increasing the macroion charge. 
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I. INTRODUCTION 

Platelike colloidal particles play a decisive, constitu- 
tive role in processes like agriculture, construction, oil 
drilling, or coating. This wide range of applicability mir- 
rors a very rich phase behavior of suspensions of platelets, 
including liquid crystalline phases, sol-gel transitions, 
and flocculation, depending on numerous material pa- 
rameters, like size, shape, or charge of the particles, as 
well as on effective, solvent mediated interactions which 
can be tuned, e.g., by the choice of the solvent, salt 
concentration, or pH-value. Whereas bulk properties of 
suspensions of charged platelets have been investigated 
for decades [3J, free interfaces between coexisting fluid 
phases in such systems have not yet been studied. Here 
we focus on this latter issue by proposing and studying 
a density-functional theory of such inhomogeneous mul- 
ticomponent systems of charged anisotropic particles. 

On the experimental side, characterization of bulk 
phases have been conducted for several model systems 
like natural clay 0,0, laponite [j H IE Jlllji EEl 
13 Elll, sterically stabilized gibbsite [IHEflOllll, 
or nickel(n)hydroxide I2CJI usin g m ethods like polar- 
ized light analysis @, TlHt Il6t ITil Il8j|. light scattering 
small-angle scattering with neutrons or x-rays 
0, [E EIIM H3 , rheological measurements 0, |H , or 
NMR jlj, [lj] . Since coexistence between bulk phases of 
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charged gibbsite platelets has already been observed ex- 
perimentally we expect that the spatially varying 
structural properties between them are also experimen- 
tally accessible using, e.g., scattering [21] or optical p2^ 
techniques. 

The theoretical description of suspensions of charged 
platelets is rather complicated due to long-ranged, 
anisotropic interactions and many different length scales. 
Interfaces and surfaces in such systems add further dif- 
ficulties induced by the partial loss of translational sym- 
metry. Under these circumstances it is advisable to start 
with simplified models. In the presented one, the parti- 
cles are modeled as hard cuboids with pointlike charges 
concentrated in their center. Furthermore, the platelet 
orientations are restricted to three mutually perpendicu- 
lar directions, which is commonly known as the Zwanzig 
model p3|. Finally, in order to gain computational ad- 
vantages, distances between particles are not measured 
by the usual Euclidean norm but by the so-called supre- 
mum norm. 

Density-functional theory [24) is a very effective 
method to investigate inhomogeneous fluid systems. It 
has recently been applied to describe suspensions of 
platelets with pure hard-core interactions near interfaces 
and surfaces [2^, l2tl |27| . Here we study platelets in- 
teracting via a hard-core plus a Coulomb potential; fur- 
ther interactions like dispersion forces are not considered, 
which corresponds to a suspension in which the indices 
of refraction between solvent and solute are matched. 
The density functional is constructed by functional in- 
tegration of the two-particle density with respect to the 
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interaction potential [24 1 . w hich is the analogue to a 
Debye-charging process |28|. The two-particle density 
can be obtained, e.g., by interaction site model calcula- 
tions [2^. l3p| or expressed in terms of the potential of 
mean force which may be approximated by effective pair 
potentials [j^ElES]- For 

reasons of computational ad- 
vantages, here we choose an extension of the even simpler 
Debye-Hiickel pair distribution function psl ] in which the 
Debye screening factor is replaced by a spatially varying 
quantity. 

In view of these approximations, this theoretical model 
can be expected to be only qualitatively correct. The 
most subtle point is the choice of the pair distribution 
function, which requires further improvements in order to 
increase the quantitative reliability. However, the other 
parts of the formalism are expected to remain valid and 
thus provide the basis for future efforts. In this sense, 
the following sections present a generic formalism for free 
interfaces in fluids of charged platelike colloids, imple- 
mented exemplarily for the above-mentioned extended 
Debye-Hiickel pair distribution function. 

The text is structured as follows. In Sec. [H] the gen- 
eral formalism with a detailed derivation of the actual 
density functional is presented. Section ITTT1 is devoted to 
bulk phase diagrams of the model. The structures of free 
interfaces between an isotropic and a nematic bulk phase 
are calculated in Sec. IIVI Section[V]discusses the current 
approach and summarizes our results. 



II. GENERAL FORMALISM 

A. Definitions 

We consider a ternary mixture of charged hard square 
cuboids with their edges required to be parallel to the 
Cartesian axes (Zwanzig model |23j ) dissolved in a di- 
electric solvent (e.g., water) with dielectric constant e. 
The solvent is treated as a continuum. For simplic- 
ity, the charges are fixed, monodisperse, and concen- 
trated in the centers of the particles. The particles of 
the first component, representing the macroions M, have 
size D M x D m x Lm, D m ^= L M , and charge Q M < 0. 
Within the Zwanzig approximation, macroions can take 
three different orientations, denoted as M x , M y , or M z 
corresponding to whether the LM-edges are parallel to 
the x-, y-, or z-axis, respectively (see Fig. QJ. The second 
component consists of salt anions A modeled as cubes of 
side length Da ■= D$ and charge Qa '■= Qs < (see Fig. 
0. Finally, the third component consists of salt cations 
C and counterions guaranteeing overall charge neutral- 
ity. They are also described by cubes with the same side 
length Dq := D$ but opposite charge Qc ■= —Qs > 
(see Fig. [Q. 

We denote as Qi(r),i e {M x , M y , M z , A, C}, the num- 
ber density at point r of the centers of macroions with ori- 
entation M Xj y iZ , anions, and cations, respectively. Note 
that the position r S V C M 3 , with V denoting the 



Qm < Qm < Q M < 




M x M y M z 



Qa = Qs<0 Q c = -Qs > 




FIG. 1: Macroions M are square cuboids of size Dm x Dm x 
Lm, Dm 7^ Lm with charge Qm, whereas anions A and cations 
C are cubes of side length Ds with charges Qs and — Qs, re- 
spectively. The pointlike charges (•) are localized in the cen- 
ters of the cuboids. The macroions can adopt three possible 
orientations M x , M y , and M z corresponding to the Z/M-edges 
being parallel to the x-, y-, and z-axis, respectively. 

system volume, is a continuous variable in contrast to 
the orientation of macroions, which varies within a dis- 
crete set. As an abbreviation, we introduce g(r) := 
(0M x (r), . . . , £>c(r)). 

The system under consideration is coupled to two 
particle reservoirs: One supplies neutralized macroions 

(chemical formula C&M, k :— — — ) and the other neutral 

Qs 

salt (chemical formula CA); /ic fc M and ficA denote the 
corresponding chemical potentials. Upon entering the 
solvent, these molecules dissociate: 

C k M — > kC Qc +M Qm , 
CA — ► C Qc + A Qa . (1) 

These equilibrium chemical reactions lead to the follow- 
ing relations between the reservoir chemical potentials 
(MCfeM and hca) and the particle chemical potentials 
(pi,i G {M x ,M y ,M z , A, C},^ Mx = A*M y = MmJ : 

Mca = + MA- (2) 

B. Density- functional theory 

The configurations of this system are characterized by 
the set of number density profiles g. The equilibrium 
states minimize the grand canonical density functional 

nm 

%] = ]T U\ Qi (v) (In ( ft (r)) - 1 - ti)+F c *{Ql (3) 
* v 
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where F ox is the free energy in excess over the ideal gas 
contribution. Here, the reduced particle chemical po- 
tentials fi* := m — In (Af ) with the thermal de Broglie 
wavelength Aj ( Am x = Am = Am z ) for particles of class i 
have been introduced. With the reduced reservoir chem- 
ical potentials 



M Cfc M : = Mc fc M- (ln(A Mxy J +fcln(A^ 
Mca ■= MCA- (ln(A|) +ln(A^)) , 
Eq. (J2J takes the form 

Mc fc M = ^Mc + Mm XiYiZ 
Mca = Mc + Ma- 



(4) 



(5) 



For given reservoir chemical potentials /ic fc M an( ^ Mca> 
the particle chemical potentials /j,* are fixed by Eq. (J5J 
and the constraint of global charge neutrality: 



d 3 r 



QiQi(r) = 0. 



(6) 



The Euler-Lagrange equations of the minimization 
problem read 



= ln(ft(r)) - Mi - Ci(r) = 



with the one-particle direct correlation function 



Ci(r) 



5F C 



(7) 



(8) 



K Q = minimizes the density functional in 

Eq. 0, the grand potential fi(T, V, Mc h Mi Mca) = 
—p(T, V, MCfcM' Mca)^ with the osmotic pressure 
p(T, V,Mc m>Mca) e Q ua l s ^[£ cq ]- Phase coexistence 
corresponds to different states with equal values of the 
pressure p, the chemical potential Mc fc M °f the reservoir 
of neutralized platelets, and the chemical potential /z£, A 
of the salt reservoir. In particular, coexistence does not 
imply equal values of the particle chemical potentials [i* , 
i £ {M x , M y , M z , A, C}. Rather, coexisting bulk phases 
give rise to a Donnan potential maintaining different 
chemical potentials [i* . 



C. Excess free energy 

The above considerations are valid for any interaction 
between the particles. In this subsection, our choice of 
the model and the resulting excess free energy F ex are 
specified. 

The interaction energy Uij(r, r') of a particle of class i 
at position r with a particle of class j at position r' com- 
prises a hard-core potential E/y(r, r'), which prevents the 
particles from overlapping, and a contribution U^Ar, r') 

due to the charges: U = U h + U c . As stated in the In- 
troduction (Sec. [U, we do not consider dispersion forces. 



The interactions between the charges are approxi- 
mated as 



QiQj 



r - r 



(9) 



where the usual Euclidean norm ||r|| 2 = \/ x 2 + y 2 + z 2 is 
replaced by the supremum norm ||r|| = max(|a;|, |y|, \z\) 
because of computational advantages. Since these two 
norms are equivalent, i.e., < ||r|| 2 < V^HrH^, 

we do not expect that the results change qualitatively 
due to this approximation. Furthermore, ||-|| -spheres 
are cubes with their edges parallel to the Cartesian axes; 
therefore, the supremum norm is the most natural and 
adapted norm in the context of a Zwanzig model for 
cuboids. 

As described in Ref. [24j , the exact relation 

,,/ " X 7 r = ^ i (r)^(r')^(r,r') (10) 



SUij (r,r 



with the pair distribution function g corresponding to the 
pair potential U can be functionally integrated along the 
path ■= U h + nU c ,n£ [0, 1], which yields 

f CX y = F^ h [g] + i£ Jd\ J dV 



y v v 
i 



ei {v) Qj {r')U^T') ItotgWM, (11) 



where F ex,h is the excess free energy corresponding to the 
pure hard-core potential U h , and denotes the (inho- 
mogeneous) pair distribution function for the pair poten- 
tial t/w. f cx > h is chosen as the fundamental measure 
functional introduced by Cuesta and Martinez-Raton 

[allil: 



F e *' h [g] := Jd 3 r $ (n(r)) 



(12) 



with the weighted densities 

n a (r) =Y,[ dV <^( r - r ')ft(r') (13) 



for a £ {0, lx, ly, lz, 2x, 2y, 2z, 3} and the excess free 
energy density 



$(ri) 



-n 



ni q n 2q n 2x n2yn2z 



ln(l-n 3 )+ > + 7- , 

f— ' - 1 - n 3 (1 - n 3 ) 



q£{x,y,z} 



2 ' 



(14) 



Due to Eq. the fy-integration in E q. (1 1 1|> may 

be interpreted as a Debye charging process [23, |37| ■ This 
motivates to approximate g^ by an expression similar to 
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the pai r distribution function of the Debye-Hiickel theory 

mtniniai: 



with functions 



fJ 



(v) 



(r,r') := cxp (-E/£(r,r')) max [0, 1 - tf§(r,r') 

x ryexp (-^^(r, r'; [g]) ||r - r']^)] . (15) 



This pair distribution function has non-negative values, 
it vanishes within the hard-core, and it approaches unity 
at infinitely large distances. Again, here the Euclidean 
norm ||-|| 2 has been replaced by the supremum norm 
IHloo- The factor y/rj in Eq. (fl5|l is introduced be- 
cause in a Debye charging process all charges Qi are re- 
placed by yfr\Qi. Furthermore, the charges Qi have to 
be interpreted as effective charges in order to reproduce 
the actual effective interactions between charged parti- 
cles within Debye-Hiickel theory |40| . 

The screening factor k in Eq. I|15l) is chosen as 

K (r,r';[^):=i(K(r;y) + «(r';[ £ ])) (16) 

where only anions and cations contribute to the screening 
113: 



«(r; [Q}) ■= ^nQl(g A (r) + g c (r)). (17) 

Our analysis rendered that, for a spatially constant 
screening factor k, the above model does not yield stable 
interfacial profiles. This led us to introduce the spatially 
varying expression in Eq. I|16|) . Alternative expressions 
for inhomogeneous screening factors are known from the 
theory of electrolytes: In Refs. [3^ . |4^| non-symmetric 
screening factors are provided whereas in Ref. |43j k is 
calculated from the mean salt density. We prefer the def- 
inition in Eqs. (|16l) and (|17f) because it is symmetric and 
the screening is determined by the salt concentration at 
the actually investigated positions. 

The expansion of £l[g = £>( hom ) +Sg] around a spatially 
homogeneous state p( hom ) in powers of perturbations Sg 
shows that the spatially homogeneous state g( h ° m > is un- 
stable with respect to spatial variations if the macroion 
charge |Qm| is sufficiently large because the second or- 
der term can become negative. Thus, the choice for 
Eq. Q15f) leads to spatially inhomogeneous bulk 
phases if the macroion charge |Qm| is larger than some 
threshold value. Here, we restrict ourselves to the case of 
spatially homogeneous bulk phases, i.e., only sufficiently 
small macroion charges are considered. 

With Eq. I|15ll . the innermost integral in Eq. I|ll(l can 
be evaluated leading to an expression 



(18) 



1 

G i:j (K, s) := - J dr] min [l, J7? (s)r?exp {-y/rjns)] (19) 



which decay for s — ► oo as 



Gij(n,s) 



dr) U[ ••(s)r?exp (-y/rjKa) 



1 

■ * J / drj 77 exp (—y/rjKs) 



QiQj 



J dC 2( 3 exp (-C«s) 



iQ j 



(20) 



where 7 denotes the incomplete gamma function Hl|45|. 
Therefore, the integrand in Eq. Hll|) vanishes at small 
distances ||r — r'Hoo due to U^, whereas it decays as 
||r — r'H^ 1 for ||r — r'|| — ► 00. In order to isolate the 
|| r — r'|| 1 asymptotics, we add and subtract unity on 
the right-hand side of Eq. (T5J, which, after insertion 
into Eq. I|ll|) . leads to the following decomposition of 
the excess free energy 



(21) 



with the electrostatic part 



\J2 J d3r / dV ftW^(r')^(r,r') 



v v 



2 J J r-r'lloo 

V V 

i / d 3 r e Q (r)^(r) 



(22) 



exp (— U£j(r, r')) (l + Gy(«;(r, r'; [g\), ||r — r'H^)) and the correlation part 
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F^M ■■= \ E / dV / dV ^) gj (r')U^(r, r') [exp (-[$.(r, r')) - 1 

+ exp(-[^.(r,rO)Gy(«(r ) r';y),||r-r'|| oo )] . (23) 



Here, the local charge density 



and the electrostatic potential 



■0(r) 



3 ' 



dV 



pQ(r') 



(24) 



(25) 



have been introduced. Note that, although the inte- 
grands in Eqs. I|22|l and l|23|) are undefined for r = r', the 
three-dimensional integrals exist due to the ||r — r'|| 



0. 



Since the integrand in 
.'i 



asymptotics for ||r — r'|| 

Eq. O decays as ||r - r'||^ 6 for ||r - r'H^ — -> oo, F™£ 
is well-defined for all finite system volumes V and the 
thermodynamic limit of the ratio 
ing to the last expression in Eq. 



F^c C orr exists. Accord- 
l|22|) . the same state- 



ments are true for the electrostatic contribution F el ' 
provided the electrostatic potential ijj is well-defined. For 
locally charge neutral systems (g Q = 0), e.g., for bulk 
phases, the latter holds because of Eq. i|25|) . In the next 
subsection it is shown that ip can also be calculated in 
systems with only lateral translational symmetry. 



where a' denotes a two-dimensional vector in the x-y 
plane. The inner integral in Eq. I|27|l leads to 



dV 



A(L) 
\x—x' I 



Ka',2- z') 



da' 8a' 1 



da'Sa'i 
\x — x'\ J a' 

\x~x'\ 

= -A\x - x'\ + 8L. (28) 

In conjunction with the global charge neutrality con- 
straint of Eq. © , Eq. (|2"T|) reduces to 



rl>(z) = -4 / dz' g Q (z')\z-z'\. 



(29) 



By differentiating twice, one finds that -0 fulfills the Pois- 
son equation: ip" = — 8g^. Furthermore, by making use 
of global charge neutrality (Eq. @J), one finds 



'•(~) = -4 / dz> Q 1(z')(z'+^ 



D. Planar geometry 

By imposing suitable boundary conditions, we consider 
only systems with translational symmetry in the lateral 
x and y directions. Hence, in the absence of spontaneous 
symmetry breaking, all densities Qi depend at most on 
the z coordinate. Since the thermodynamic limit of glob- 
ally charge neutral systems of Coulomb interacting hard 
particles exists |46j, i.e., the bulk free energy density de- 
pends neither on the shape nor on the boundaries of the 
system volume V, the following system volumes of size 
2L x 2L x L are considered in the limit L — > oo: 



V{L) :=A(L) x 



L L 
2"' 2 



(26) 



where A(L) is a square in the x-y plane of side length 2L 
with periodic boundary conditions. 

The electrostatic potential in Eq. (|25|l can be ex- 
pressed as 



i>{z) = / dz' g Q (z') / dV 



A(L) 



||(a',z-z')ll 



-, (27) 



4 / dz' Q %')(±-z' 



(30) 



Thus, the density functional in Eq. takes the final 
form 

L_ 

2 

n[g] = 4L 2 E J ^ Qi (z) (ln( Qi (z)) - 1 - M | 



+ gQ^C*) 1 + F + F%°\£ (31) 

which has to be minimized under the constraint of global 
charge neutrality (Eq. ©)• This is carried out by nu- 
merically solving the Euler-Lagrange equations (Eq. Q) 
with a Picard-iteration scheme on a one-dimensional grid. 
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In the isotropic and nematic bulk fluid, the densities g 
are spatially constant. In this case, the Euler-Lagrange 
equations Eq. Q comprise five coupled equations: 

, , fbulkK *(bulk) h,(bulk) c(bulk) _ „ r „„. 
m {t>i ) Mi H c coir,i ~ u - V ^) 

where the chemical potentials ^*( bulk ) fulfill Eq. JSJ and 
the local charge neutrality condition Qig^^ = 0). 
The electrostatic contribution — Cgj^ bulk ^ = Qi^ bulk ' is 
absent in Eq. (|32|) because the electrostatic potential ip 
vanishes in locally charge neutral systems (see Eq. I|25|l). 

For determining the number density profiles at free in- 
terfaces between coexisting bulk phases Bi and B2, the 
Euler-Lagrange equations 

In (g t (z)) - M * + Q^(z) - $(z) - c c co „Az) = (33) 

are to be solved with the boundary conditions 



L \ (81 
Qi\ z = ) = Qi 



'->> \ z = 2 I = L ' 



(B2) 



(34) 



In order that for z = — ^ and z = U Eq. (|33|) reduces 
to Eq. (|32|) for B\ and B2, respectively, one has the 
requirements 



M* +Q i i> 



L 



*(Bi) 



-m* + Qi4> (§) = -^ B2 

Using Eq. I)3()[l. one readily concludes 



(35) 



(36) 



and 



_ (37) 
%1>d is known as Donnan potential |47| between the two 
bulk phases B\ and Bi- It maintains a density gradient 
of the mobile particles at the interface between two coex- 
isting bulk phases. Its definition given above is unique, 
i.e., the rightmost expression is in fact independent of i 
due to Eq. JSJ), e.g., 



.(Si) _ *(«0 

Ma Ma 



(^A"Mc (6l) )-(^CA-Mc (B2; 



-Qc 



*(Bi) *<B 2 

-Mo 



(38) 



III. BULK FLUID 



As a first step in the investigation of the density func- 
tional developed in the last section, bulk phase diagrams 



are determined for various macroion charges Qm by solv- 
ing the bulk Euler-Lagrange equations 

The macroion and the salt number densities are given 
by Qm '•= QMx+QMy + QM* and g s := g A , respectively. In 
order to detect the formation of liquid crystalline phases 
of the macroions, the equilibrium nematic order param- 
eter for the director oriented relative to the z-direction, 



s M 



3 £>M Z 
2 £>M 



-,1 



and the equilibrium biaxial order parameter, 



Qm 



QM 



(39) 



(40) 



have been determined. The definition of sm agrees 
with the well-known scalar liquid-crystal order param- 
eter S — (P 2 (cos^)) = |((cos$) 2 ) — \ because within 
the Zwanzig model only macroion orientations M z par- 
allel (cost? = 1) and M XjY perpendicular (cosi? = 0) 
to the z-axis are possible, sm vanishes in an isotropic 
phase (gM x = gm y = f?M z ), whereas it is positive in 
a nematic phase with director parallel to the z-axis 
(pm z > £>M x ,£>M y )- A discrimination of the orientation 
M z leads to negative values of sm- 

It turned out that the biaxial order parameter qyi van- 
ishes throughout the whole inspected range of reduced 
chemical potentials ^c k M an d A^caj whereas the nematic 
order parameter sm indicates either an isotropic fluid 
(sm = 0) or a nematic fluid (% > 0). 

Figure El shows phase diagrams for the parameters 
(compare Fig. UJ) D M = 20 £ B » 14 nm, L M = 
l B ~ 0.72 nm, D s = t B » 0.72 nm, Q s = -e with 
Qm = 0, Q M = 0.25Q S , Q M = 0.5Q S , Qm = 0.75Q S , 
and Qm = Qs m terms of the macroion packing fraction 
T]M = £>m-D m L m and the salt density gs. 

One isotropic phase (I) and one nematic phase (N) 
are found separated by first-order phase transitions. 
Whereas for coexisting phases r]M is always smaller in the 
isotropic phase than in the nematic phase, gs of coexist- 
ing phases is higher in the isotropic and lower in the ne- 
matic phase. A similar displacement of salt from regions 
of large concentrations of charged macroions is known as 
Donnan effect 01- Whereas the original Donnan effect 
has been discovered in systems subdivided by membranes 
which arc impermeable for macroions, here the density 
difference of the macroions occurs due to two coexisting 
bulk phases. As for the case of membrane equilibrium, 
here a Donnan potential 4>d (Eq. (|S3) maintains the 
density gradients between the coexisting phases. Figure 
13 shows its dependence on the macroion charge Qm and 
the salt density gs- ipD decreases with increasing salt 
density gs- This tendency is intuitively expected as the 
Donnan effect becomes more pronounced with increasing 
macroion charge whereas increasing the salt density gives 
rise to a stronger screening of the macroion charge. For 
fixed salt density gs well below w 0.1 mM, ipu decreases 
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FIG. 2: Bulk phase diagrams of mixtures of plate-like macroions (Dm = 20 Ib, Lm = Ib) and monovalent salt (Ds = is, 
Qs = -e) for macroion charges Q M = (a), Q M = 0.25Qs (b), Qm = 0.5Qs (c), Qm = 0.75Qs (d), and Qm = Qs (e) in 
terms of the platelet packing fraction t]m = ^m^m^m an d the salt density gs |34| . Note that gs = ImM for Ds = 0.72 nm 
corresponds to a salt packing fraction QsDg = 2.2 ■ 10 -4 . Coexisting states are connected by (non-horizontal) tie lines. The 
model exhibits one isotropic (I) and one nematic (N) phase separated by first-order phase transitions. The salt density of 
coexisting phases is higher in the isotropic and lower in the nematic phase (Donnan effect) as can be inferred from the negative 
slope of the tie lines. For increasing macroion charge, the isotropic and the nematic binodals are shifted to larger macroion 
packing fractions i]m- With increasing salt density gs, the isotropic-nematic binodals for systems of charged macroions bend 
towards lower values of the macroion packing fraction. 



with increasing |Qm|, whereas this behavior is reversed 
for fixed salt density gs above « 0.1 mM. 

Upon an increase of the macroion charge |Qm|, the 
isotropic and the nematic binodals are shifted to larger 
values of the macroion packing fraction t/m- This may 
be qualitatively understood by introducing the notion 
of an effective shape, which, in the present case, for 
macroions is given by a hard core surrounded by a soft 
Ij-H^-sphere, i.e., a cube, with its linear extension pro- 
portional to Q M due to the pairwise Coulomb repulsion. 
For small macroion charges, the effective shape is still 
platelike whereas for highly charged colloids, the effec- 



tive shape tends towards a cube leading to a shift of the 
two-phase region to larger macroion packing fractions. 

For fixed macroion charge as well as particle shape and 
increasing salt density gs , the isotropic-nematic binodals 
in Fig. [2] bend towards smaller macroion packing frac- 
tions. This behavior is expected intuitively, because high 
ionic strength causes strong screening which in turn leads 
to effectively quasi-hard platelets (see Fig. EJa)). 
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FIG. 3: Donnan potential ipn (see Eq. (1371 "1 between the 
coexisting nematic (N) and isotropic (/) bulk phases in mix- 
tures of platelike macroions and salt (see Fig. |2J. For fixed 
macroion charge Qm, the Donnan potential decreases with 
increasing salt density gs- For salt densities gs ~ 0.01 mM, 
ipD decreases with increasing |Qm|, whereas for salt densities 
gs ~ lmM, ipo increases with increasing |Qm|. 



IV. FREE INTERFACES 

Based on the bulk properties provided in the previous 
section, we are now able to calculate the density profiles 
at the free interfaces between the coexisting isotropic and 
nematic phases by solving the spatially varying Euler- 
Lagrange equations (Eq. The density and order 

parameter profiles corresponding to the parameters used 
in the previous section (see also Figs. |21 and |3J) and to 

nematic bulk salt density gf ] = 2.2 • 10~ 5 = 0.1 mM 
[34j are depicted in Figs. |t| and respectively. The 
interface position z = is chosen such that £>m(0) = 



(N) 



First, and most important, the formalism described in 
Sec. [n] renders stable free interfaces between coexisting 
bulk phases. This should be regarded as an accomplish- 
ment which can be traced back to using the spatially 
varying screening factor k introduced in Eq. (|16(l : trials 
with spatially constant n were not successful. 

For a given density profile qm(z) (Fig. 0J), the corre- 
sponding interface width £ is defined as the spatial dis- 
tance between the loci, where the tangent at the density 
profile at position z — reaches the values of the nematic 

bulk density q^P and the isotropic bulk density g M \ re- 
spectively. This interface width £ decreases monotoni- 
cally with increasing macroion charge from £ = 1.8Z?m 
for Qm = 0.25Qs to £ = 1-3Dm for Qm = Qs (see 
Fig. 2J. The Debye length k _1 decreases monotonically 
from = 1.5£>m for Qm = 0.25Qs to = Dm 
for Qm = Qs- Finally, the bulk correlation lengths £ of 




FIG. 4: Macroion density profiles at the free interface be- 
tween coexisting isotropic (/) and nematic (iV) phases for the 
mixtures of platelike macroions and salt as studied in Fig. 
II at nematic bulk salt density g { s N) = 2.2 ■ 10 -5 = O.lmM 
|34| . The interface position z = is chosen such that 

£>m(0) = \ ^f?M + • The interface width £, inferred from 
the slope of the density profiles at the interface position z — 
(see main text), the Debye length ft -1 , and the bulk correla- 
tion lengths £ , inferred from the exponential decay of the den- 
sity profiles, decrease monotonically with increasing macroion 
charge \Qm\- The same behavior is found for the interfacial 
tension 7. See also Fig. 



the coexisting isotropic and nematic bulk phases, inferred 
from the exponential decay lengths of qm(z) — £? M , also 
decrease monotonically upon increasing |Qm| and the val- 
ues are by and large equal to those of K . 

The nematic order parameter profiles sm(z) (Fig. |3J) 
interpolate almost monotonically between sm(oo) > 
in the nematic bulk phase (./V) and sm(oo) = in the 
isotropic bulk phase (/). Note that sm has been de- 
fined for a director in z-direction, i.e., platelets on the 
nematic side (z < 0) are preferably oriented parallel to 
the free interface. At a fixed position on the nematic side 
{z < 0), sm(z) decreases with increasing macroion charge 
|Qm|- This behavior is consistent with the picture of an 
increasingly isotropic effective shape, introduced in the 
previous section. 

The charge density profiles sr{z) (Eq. (J23J) displayed 
in Fig. show deviations from local charge neutral- 
ity within the interfacial region — 4Z?m ^ z < 4£>m- 
A negative charge density occurs on the nematic side 
(N) and a positive charge density on the isotropic side 
(I). Such a local charging is necessary for the appear- 
ance of the non- vanishing Donnan potentials ipo shown 
in Fig. [3J The full electrostatic potential profiles ip(z) 
are depicted in Fig. \7\ They increase monotonically from 
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FIG. 5: Macroion nematic order parameter profiles sm (Eq. 
1391 ) at the free interface between coexisting isotropic (7) 
and nematic (iV) phases for mixtures of platelike macroions 
and salt as shown in Fig. [5] at nematic bulk salt density 
g ( s N) = 2.2 ■ 1(T 5 = 0.1 mM The position z = is fixed 
by the choice qm{0) = | + £>m^ (see Fig. 0J. On the 

nematic side of the free interface, most of the macroions lie 
parallel to the interface. For z < 0, sm(z) decreases with 
increasing macroion charge |Qm|. The inset shows a detailed 
view of the steepest portions of the nematic order parameter 
profiles in the range sm(z) € [0.2, 0.3] indicated by the frame. 



the macroion-rich nematic phase N to the macroion-poor 
isotropic phase /, maintaining the density gradients oc- 
curring in the interface region. The potential difference 
■0(oo) — ip(—oc) equals the Donnan potential ipD (see Eq. 
I|37t and Fig. E| . 

The interfacial tensions 7 of the interfaces shown 
in Fig. 21 decrease monotonically from 7 = 1.33 • 
10~ 5 = lOSnN-m- 1 for Q M = 0.25Q S to 7 = 6 • 
10~ 7 = 5nN • m _1 for Qm = Qs 0, which are compara- 
ble to experimental findings for laponite suspensions |16| . 
The corresponding wetting parameters u> = (47T7^ 2 ) -1 
[48j are in the range 6 . . . 340. If these values for 7 are in- 
deed so small, the free isotropic-nematic interfaces are 
expected to be strongly affected by capillary wavelike 
fluctuations which are not captured by the present the- 
ory. 

Calculating bulk phase diagrams and density profiles 
for macroion charges Qm > 1.25Qm along the lines de- 
scribed above lead to unphysical results like, e.g., nega- 
tive interfacial tensions. The reason for this phenomenon 
is that the bulk phases determined in Sec. 1 1 1 1 1 are as- 
sumed to be spatially homogeneous, whereas it can be 
shown that the equilibrium bulk states are spatially in- 
homogeneous for sufficiently large macroion charges (see 
Subsec. HTT1 . 



FIG. 6: Local charge density profiles (Eq. ((2U) at the free 
isotropic-nematic interfaces shown in Figs. 21 and lousing the 
same parameters and line code. Whereas global charge neu- 
trality holds, deviations from local charge neutrality occur 
near the interface with a negative charge density on the ne- 
matic side (JV) and a positive charge density on the isotropic 
side (I). The corresponding electrostatic potential profile is 
shown in Fig. 

V. DISCUSSION AND SUMMARY 

The numerical results for the bulk systems and the 
free interfaces presented in Sees. IHII and HVI respectively, 
are in good qualitative agreement with intuitive expec- 
tations. Within the density functional theory described 
in Sec. [H]one gains access not only to interfacial density 
profiles at free interfaces between coexisting bulk phases 
but also to local charge densities and electrostatic poten- 
tial profiles. 

Since here we have been interested in spatially homoge- 
neous (isotropic or nematic) bulk phases, only very small 
macroion charges (|Qm| < e) have been considered. De- 
termining phase coexistence for larger platelet charges is 
computationally more demanding because spatially inho- 
mogeneous bulk phases are involved. Note that, accord- 
ing to Subsec. Ill CI charges have to be interpreted as 
effective charges. 

A further difficulty related to the choice of the pair 
distribution functions g£- used in Eq. (jll|l may appear 
if, as for Eq. Q15[). the effective shape of the macroions 
becomes more and more isotropic upon increasing the 
macroion charge: In this case, the two-phase region be- 
tween isotropic and anisotropic phase is shifted to un- 
rcalistically large packing fractions. Therefore, the pair 
distribution functions must be chosen properly in order 
to yield anisotropic effective macroion shapes up to large 
macroion charges. Unfortunately, deriving analytical ex- 
pressions for pair distribution functions of platelike par- 
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FIG. 7: Electrostatic potential profiles at the free isotropic- 
nematic interfaces shown in Figs. 14161 The potential differ- 
ence between the nematic and the isotropic bulk is given by 
the Donnan potential tpn displayed in Fig. [3] it is indicated 
for Qm = Qs- 



tides with inhomogeneous charge distributions is still a 
big challenge. 

In conclusion, the density- functional theory of charged 
platelike particles developed here can be regarded as a 
first step to understand qualitatively free interfaces be- 



tween isotropic and nematic bulk phases of suspensions of 
charged platelets and salt. Here, we have constructed a 
density functional for charged platelike particles and ap- 
plied it to a ternary mixture of platelike macroions and 
salt ions (Fig. 0) in the bulk and at free interfaces be- 
tween coexisting isotropic and nematic phases. For suffi- 
ciently small macroion charges, the bulk phase diagrams 
in terms of densities (Fig. |2J) exhibit one isotropic phase 
and one nematic phase. For increasing macroion charge 
and fixed salt densities, the two-phase coexistence region 
is shifted to larger macroion packing fractions. For fixed 
macroion charge and increasing salt density, the limit of 
quasi-hard platelets is approached. The Donnan poten- 
tial between coexisting phases (Fig. can be expressed 
in terms of the particle chemical potentials gained from 
bulk structure calculations. Density and nematic order 
parameter profiles at free interfaces between isotropic 
and nematic phases at coexistence (Figs. 0]and|5Jl show 
non-monotonic behavior. The value of the nematic or- 
der parameter in the nematic bulk phase decreases upon 
increasing the macroion charge. The width of the inter- 
face and the bulk correlation lengths are approximately 
given by the Debye length. The interfacial tension de- 
creases upon increasing the macroion charges. Electri- 
cally charged layers form at the free interface (Fig. [fJJ). 
The corresponding electrostatic potential profiles (Fig. 
0) exhibit monotonic behavior. Investigations of spa- 
tially inhomogeneous bulk phases are necessary in order 
to apply the theory to larger macroion charges. Improve- 
ments of the present theory call for more accurate analyt- 
ical expressions for the pair distribution function between 
charged platelike particles. 
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